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We derive an explicit expression for the coupling constants of individual eigenstates of a closed 
billiard which is opened by attaching a waveguide. The Wigner time delay and the resonance 
positions resulting from the coupling constants are compared to an exact numerical calculation. 
Deviations can be attributed to evanescent modes in the waveguide and to the finite number of 
eigenstates taken into account. The influence of the shape of the billiard and of the boundary 
conditions at the mouth of the waveguide are also discussed. Finally we show that the mean value 
of the dimensionless coupling constants tends to the critical value when the eigenstates of the billiard 
follow random-matrix theory. 
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I. INTRODUCTION 
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During the last years quantum chaotic scattering was 
a field of intense research. A great deal of the results ob- 
tained was based on the projection operator formalism 
due to Feshbach, Weidenmiiller and others [[jj g|. In 
this approach the scattering system is decomposed into 
a closed subsystem described by the internal Hamilto- 
nian H ln with discrete bound states n = 1 . . . N and a 
continuum of external scattering states labeled by the 
energy E and an index A = 1 . . -A.(E) corresponding to 
different open scattering channels. The coupling between 
the internal and external subsystems is then incorporated 
by an operator with matrix elements W nt \{E). The S- 
matrix of the complete system can be expressed in terms 
of these matrix elements and the Hamiltonian H m . This 
relation can be cast into the form 
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Here, S and K are energy-dependent square matrices of 
dimension A x A, and W has dimension N x A. While 
this setting is very general, the tools developed for the 
subsequent analysis of the properties of the S-matrix re- 
quire additional assumptions. In particular, the energy 
thresholds for the opening of new scattering channels are 
usually neglected. As a consequence the energy depen- 
dence of the coupling matrix W(E) can be considered 
weak and then the S-matrix ([!]) can be rewritten in terms 
of an N x N effective non hermitian Hamiltonian H e g 
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This canonical formalism |y, g, |[ for expressing the S- 
matrix is sometimes denoted as the Heidelberg descrip- 
tion of scattering, and we use this name to distinguish it 
from an S-Matrix obtained directly, i. e. without refer- 
ence to any auxiliary closed system. 

When the internal Hamiltonian in (Eh describes a 
chaotic system, it is justified to replace it by a random 
matrix |4|], and by performing an average over the ap- 
propriate ensemble a statistical theory for the S-matrix 
is obtained which allows to calculate quantities of inter- 
est such as correlation functions or the distribution of 
Wigner delay times and resonance poles || |[ @, ||. It 
was found that the results of such an approach are to 
a large extent independent of the detailed structure of 
the matrix W, but they do depend on the dimensionless 
mean coupling strength 
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Here, D is the mean energy level spacing of the inter- 
nal subsystem, and the average (. . .) n .\ is taken over the 
internal states n and all open scattering channels A. 

For example, when the coupling constant (|5|) exceeds 
the critical value g = 1 and the number of scattering 
channels is small compared to the total number of states, 
a counterintuitive shrinking of the widths of most reso- 
nances with increasing coupling is observed || [l0| pd| . 
For each attached scattering channel only one of the res- 
onance widths grows further with the coupling g. The 
resulting redistribution of S-matrix poles was coined res- 
onance trapping. 

Chaotic billiards with attached waveguides are con- 
sidered as paradigm for chaotic scattering They 
are relevant as theoretical models for understanding the 
transport properties of mesoscopic semiconductor struc- 
tures p3| or experimental results on microwave scattering 
in flat resonators |L4|, 15 . Also signatures of resonance 



2 



trapping were recently observed in billiards both, numeri- 
cally jl(| [I?]. [Is) and in microwave resonator experiments 

©■ 

However, to our knowledge there is no theory which 
maps a given billiard to an effective Hamiltonian with 
overcritical coupling, thus really establishing a connec- 
tion between the numerically and experimentally ob- 
served phenomena and the results on resonance trapping 
obtained within the formalism ([!])- (|5|). 

Motivated by this situation, it is the purpose of the 
present paper to discuss the application of the Heidel- 
berg approach to open billiards in some detail and to 
answer questions such as: How can (|l|), (||) be derived 
for a billiard, what kind of approximations are involved 
and what is the resulting expression for the coupling con- 
stants W n<x (E)? What is the influence of the choice of 
the internal subsystem which is not unique for a given 
scattering system? 

Using the expression for the coupling constants W n< \ 
to be derived in section II we will then address the effec- 
tive coupling constant for a typical chaotic billiard with 
an attached waveguide. We show in section III that in 
the semiclassical regime, and when no tunneling barriers 
obstruct the waveguides, the coupling strength is fixed at 
the critical value g = 1, independently of the size or the 
precise geometry of the billiard and of its openings. A 
numerical verification of our results is contained in sec- 
tion IV, followed by a short discussion of the implications 
of our findings. 

II. COUPLING CONSTANTS FOR 
INDIVIDUAL LEVELS 




FIG. 1: A scattering system consisting of an infinite waveg- 
uide and a cavity is shown with bold lines. Various possibili- 
ties to add a wall and obtain a closed billiard are shown with 
thin lines (solid and dotted). 

We consider a situation as shown in Fig. [j]. A scat- 
tering system is formed in two dimensions by an infinite 
waveguide of width b and and an arbitrary cavity. In- 
side the system the potential is identically 0. We set 
h = 2m = 1 and E = k 2 such that the stationary 
Schrodinger equation reduces to the Hclmholtz equation 

(A + k 2 )*(x,y) = 0, (6) 

with A = d 2 /d 2 +d 2 /d 2 . On the boundary of the scatter- 
ing system (bold solid line in Fig. flh we require Dirichlet 



b. c. W = 0. This boundary condition and also the 
precise geometry of the system are by no means essen- 
tial, the following generalizes e. g. immediately to a cav- 
ity with more than one attached lead or Neumann b. c. 
d/d n ^ = 0. In Fig. |l| we have shown several possibilities 
to define a closed billiard which corresponds to the scat- 
tering system in question (solid and dotted thin lines). 
We will restrict the discussion to the case shown with a 
solid line: We require that the boundary of the internal 
system is located inside the attached waveguide and that 
it consists of a transversal straight line on which cither 
Dirichlet or Neumann boundary conditions are imposed 
to close the system. Clearly, even this restriction makes 
the correspondence between the scattering system and 
the auxiliary internal system not unique because the ex- 
act position of the closure along the waveguide is variable. 
We use coordinates where this closure is at x = while 
the matching between waveguide and cavity is at x = —L 
(L > 0). 

In the region of the attached waveguide (x > —L) we 
can decompose any function into transversal modes 
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because these functions form a complete and orthonormal 
basis on the interval (0, b) according to 

oo 

E^Wa(!/')=%-!/') (0<y,y'<b) (8) 

A=l 

and 

[ dy<t> x (y)<t> xl (y) = 6 xx ,. (9) 
Jo 

The most general solution of the Helmholtz equation is 
a superposition of the scattering states ^\(x,y) which 
consist of a single incoming wave in transversal mode A 
and the corresponding outgoing modes given by the S- 
matrix of the system 

e -ik>,x +ik x ,x 

v x (x, y ) = Mv)—t^ + E 5 *<A Mv) -77— ■ (10) 

V^A y V^A' 



The longitudinal wave number k\ = \/fc 2 — (\ir/b) 2 is 
real for A < A = [kb/ir], where [. . .] denotes the integer 
part. These modes are called open or travelling, and the 
A x A matrix S\>\ corresponding to the open modes is the 
unitary S-matrix we are interested in. For A > A the mo- 
mentum along the waveguide is imaginary. These modes 
are called closed or evanescent. In the scattering state 
with A < A the evanescent outgoing modes describe expo- 
nentially decaying contributions which modify the wave 
function in the vicinity of the mouth of the waveguide. 
Evanescent incoming modes are exponentially increasing 
into the waveguide and thus unphysical for the scattering 
system. 
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When the evanescent modes are included, the S-matrix 
becomes an infinite-dimensional operator which is no 
longer unitary. It is possible to construct an eigenstate 
of the closed billiard by a superposition of the scattering 
states ( |l0| ) including evanescent modes: Suppose S(E) 
has an eigenvalue unity at some energy E n — k^ and let 
0"n,x be the components of the corresponding eigenvector. 
Then the linear combination of scattering states 

„( N ) 



7 (N) 
L n,X 



4>x(y) COs(fc n ,A x) (11) 



{—L < x < 0) satisfies Neumann b. c. at x = and it 
is thus indeed an eigenfunction of the billiard because it 
satisfies (||) and the remaining boundary conditions by 
construction. The normalization of the S-matrix eigen- 
vector a n x in ( fill) is not unity but rather determined by 

the normalization of the billiard eigenfunction ^^(a;, y). 
If Dirichlet boundary conditions are required at x = the 
same argument can be repeated for an eigenvalue —1 of 
the S-matrix and we have 

„(°) 
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(— L < x < 0). Consequently the spectrum of the billiard 
closed with Neumann or Dirichlet b. c. can be found from 
the secular equation 



det(lTS(E)) = 0, 



(13) 



which was first derived by Doron and Smilansky p9[ . 
In a sense we will in the following invert this so-called 
scattering approach to the quantization of billiards: we 
will express the S-matrix in terms of the eigenvalues and 
eigenf unctions of the closed system. 

For this purpose consider the Green function of the 
closed billiard which is defined as the resolvent of — A in 
the space of functions which satisfy the boundary condi- 
tions of the billiard. In position representation this defi- 
nition can be expressed by the inhomogeneous Helmholtz 
equation 



(A + fc 2 )G(r;r',fc) = S(r - r') . 



(14) 



In the eigenbasis of the billiard the Green function reads 

" *:(r')*n(r) 



G(r;r',fc) = ^ 



L.2 _ U2 



(15) 



which can be verified using A*f? n (r) = — k^ n (r) and the 
completeness of the functions ^f n inside the billiard. For 



r and r' inside the waveguide we can expand the Green 
function with respect to the transversal modes 4>\(y) and 
find as the general form of a solution of (|l4|) 



G(r,r',fc) = jJ2 

XX' 



4>x(y) fix'iy') 
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Indeed, the first term inside the brackets gives rise to one 
particular solution of the inhomogeneous equation (H4) , 
while the second term with the unknown matrices G^~, 
G H , G h and G represents the most general solution 
of the corresponding homogeneous Helmholtz equation 
(^). The unknown coefficients must be determined such 
that the Green function satisfies also the boundary con- 
ditions inside the cavity and on the transversal closure 
of the waveguide. For this purpose assume first x > x' 
and consider r as a fixed parameter. Then G(r; r', k) as a 
function of r' should satisfy the homogeneous Helmholtz 
equation with the boundary conditions of the scattering 
system, i. e. it can be written as a superposition of the 
scattering states VP x(r') defined in Eq. (|10|). On the other 
hand, when r' is fixed, the Green function as a function of 
r satisfies the boundary conditions (Neumann or Dirich- 
let) at x — where the closed billiard is separated from 
the waveguide by the additional straight wall. Thus, it 
must be a superposition of the functions 



^N/D (r) = ^y) (eifeAX±e _ ifeAX)! 



(17) 



which are in fact the scattering states for a semi-infinite 
waveguide with Neumann or Dirichlet b.c. at one end. 
Consequently, the Green function has the form 



G(r; r', k) = ± £ < D (r) gw (k) * A , (r') : 



(18) 



A, A' 



with another set of undetermined coefficients g\\>. Ex- 
panding (|l8|) into transversal modes and comparing to 
( filf ) we obtain 



C++ = gS G+- =g-I 
G-+ = ±gS G— = ±g 



(19) 



We can now repeat this argumentation under the op- 
posite assumption x < x' and find again the relations 
(pi) but with G+- and G~+ exchanged. This can be 
regarded as a consequence of the symmetry of the Green 
function with respect to its two arguments which in 
turn follows from time-reversal symmetry. We conclude 
G^ = G ^ = g — I = ±<7 S and hence 

g(k) = (I T S(k))- 1 . (20) 

Note that g(k) and thus the Green function diverges as 
expected at the solutions of the secular equation (|l3|), 
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i. e., when k corresponds to an eigenvalue of the closed 
billiard. 

Using Eq. ( p0| ) we can now directly relate the S-matrix 
to the transversal expansion coefficients of the Green 
function at the closure of the billiard. We define for x, x' 
inside the waveguide 



G\,y(x,x') = / dydy / (j)x(y)G(x,y;x',y / ;k)(l>x'(y') 
Jo 

(21) 

and 
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\/k\k\i d x d x 
and derive from (O) using (|2 



I + S 

iK W = (Neumann b. c. at x = 0) (24) 

1 O 



ii\T (D) = J — ^ (Dirichlet b. c. at x = 0) . (25) 

Obviously, K ( N ) = for Dirichlet b. c. and K^) = Q for 
Neumann b. c. at x = 0. Eqs. (pi|), ( f25| ) can be inverted 
and yield 



5' 



/_iif(N) /-iif(D) 



J + ii^W / + iX( D )' 



(26) 



which is now in the form of (|T|) (for Neumann b. c. up 
to an irrelevant constant phase). We can now proceed to 

determine the corresponding coupling constants W^/ T> ' > 
by representing the K-matrix in the eigenbasis of the bil- 
liard. From (0) we have 



K X \> = TT 2^ -p Jo~ 



(27) 



n=l 



and from @, @, (||) we find 
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where we have introduced the projections 



*1 N a /D) (*) = / ^ (y) *i N/D) (x, y) (30) 



of the eigenfunctions of the closed billiard onto the 
transversal modes of the waveguide. The values of the 



coup ling constants follow from comparing 
and (||) 
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(32) 



This form of the dependence of the coupling constants on 
the internal wave functions is not surprising: also within 
perturbation theory the coupling depends on the value of 
the wave function at the point where the system is opened 
or on its normal derivative for Neumann and Dirichlet 
boundary conditions, respectively. However, in the situ- 
ation we consider perturbation theory is not applicable 
and, in particular, the precise value of the prefactor in 
the coupling constants (pi|), ( j32|) could only be obtained 
from the derivation given in this section. 

The representation (|2q ) of the S-matrix in terms of 
the K-matrices (psj), (|29|) is exact, when all transversal 
modes are included. However, usually one is interested 
only in the A x A unitary part of the S-matrix, and this 
is only approximately given by (|26|) when the K-matrix 
is restricted to open modes. In [|19| the effect of this 
so-called semiquantal approximation for the accuracy of 
eigenvalues within the scattering approach to quantiza- 
tion was investigated numerically. It becomes negligible 
when the energy E is sufficiently far from the threshold 
for the opening of a new channel. Under this restriction 
we can consider as the coupling constants cor- 

responding to the unitary part of the S-matrix. When 
E approaches a threshold, Eq. (|J) breaks down, since 
the energy-dependence of the coupling constants can no 
longer be neglected. We will not consider this case here. 



III. THE MEAN COUPLING STRENGTH 

Given the explicit values (j3l|), (|3^ ) for the coupling 
constants between individual states and individual scat- 
tering channels we are now going to derive an estimate for 
the dimensionless coupling strength (^) in the semiclas- 
sical limit and neglecting evanescent modes. Since the 
concept of a mean coupling strength is not well-defined 
for infinitely many internal states, the internal Hamilto- 
nian H m entering (0), (|J) should for this purpose be cut 
to some finite matrix including only states which are close 
enough in energy k n ~ k. In particular this means that 
we can replace the momenta along the waveguide k n ^\ 
in Jil]), ( [l2] ) by their on-shell values k\. The resulting 
approximate expansions of the billiard eigenfunctions are 
projected onto the transversal modes according to ( |30| ) 
and inserted into (|3l|), (|32|) which simplify to 



W. 
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(33) 
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At this point it is necessary to determine the average 



d=l 



magnitude of the coefficients a 



(N/D) 



We assume that 



the classical dynamics of the billiard is chaotic. In the 
semiclassical limit this means that the quantum ergodic- 
ity theorem applies to the eigenstates of the billiard, i. e. 
in particular the probability density integrated over an 
arbitrary region of the billiard tends to the relative area 
of that region. Applied to the part of the billiard inside 
the waveguide we find 



bL 
T 



dy f dz|<4 N )(z,y)| 2 

J —L 
1-0 |„ |2 

- cos 2 (k\x) 

-L 



1 



£ rkb/w / 
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where A denotes the total area of billiard. In the sec- 
ond line we have inserted the normal mode decomposi- 
tion into \$F^ (x, y)\ 2 . The orthonormalization of the 
transversal modes was then used to restrict the result- 
ing double sum over modes to diagonal terms. In the 
third line (cos 2 ) = 1/2 was used ((sin 2 ) = 1/2 in the 
completely analogous calculation for Dirichlet b.c), and 
the sum over modes was approximated by a continuous 
integral. This is justified when the number of modes is 
large, i. e. in the semiclassical limit. Using the resulting 
constraint on the normalization of the coefficients of the 
billiard eigenfunctions in the transversal basis implied by 

@,(|a2 D) | 2 )=4/Awefind 



W, 



(N/D) |2 



4 

Att 



(35) 



According to (||) the average coupling between the in- 
ternal states and the continuum must be normalized by 
the mean level spacing D of the billiard which is the only 
independent energy scale of the system. To leading semi- 
classical order we have Weyl's law D = Ait/ A |21| which 
finally results in 



9 = 1. 



IV. NUMERICAL RESULTS 



(36) 



To check the validity of the Heidelberg approach for a 
quantum billiard we have performed direct numerical cal- 
culations for a Sinai billiard connected to a single waveg- 
uide (see Fig. ||). First, we evaluated numerically the 
1300 lowest eigenvalues and eigenvectors of the closed 
system with Dirichlet and Neumann b.c. at the bound- 
ary segment to which the waveguide was attached. Us- 
ing the expressions (|30|)-(|3^) we calculated the elements 




L=12 



FIG. 2: Geometry of the billiard. Four different possibilities 
to attach the waveguide are used. 



of the coupling matrix W^J D ^ . The first observation is 
that the distribution of the numerically obtained values 
of wj^/^ is fairly close to a Poisson distribution, that 
is expected for a fully random coupling. Hence the el- 
ements of the coupling matrix W evaluated with help 
of the formulae (30)-(p2|) resemble the random coupling 
standardly used in the Heidelberg approach. 

Knowing W and using Eqs. (|l|), (0) we evaluated — 
as the next step — the S'-matrix and compared it to the 
S-matrix obtained by a direct method based on the nu- 
merical solution of the underlying Schrodinger equation 
(see Rcf. ]22j for details). In order to visualize the dif- 
ferences between these two S-matrices we compared first 
the corresponding Wigner-Smith time delays 



(37) 



where M is the number of the open channels inside the 
waveguide. The results are plotted in Fig. |[ 

The figure shows an outstanding agreement of the two 
approaches for low energies. The Heidelberg approach 
describes the positions as well as the widths of the nar- 
row resonances with high accuracy. The agreement is 
good even close to the threshold energies of the individ- 
ual channels. For higher energies, however, the difference 
between the two time delay functions increases due to 
the limited number of internal states included into the 
evaluation of the Heidelberg S'-matrix. Similar results 
(not displayed) were obtained also with the waveguide 
attached to the boundaries No. 2, 3 and 4 respectively. 
In all these cases the Neumann matching procedure was 
used. 

On the other hand, for Dirichlet matching (Eq. [32]) the 
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FIG. 3: The Wigner- Smith time delay obtained for the Hei- 
delberg S-matrix with Neumann b.c. (dashed line) compared 
with the result of a direct evaluation (solid line). In the dis- 
played case the waveguide was attached to the boundary No. 
1. 



results change drastically. In this case the agreement is 
not good even for low energies. This may seem surpris- 
ing because the derivations of the previous sections were 
entirely parallel for Neumann and Dirichlet b. c. How- 
ever, an important difference is hidden in the convergence 
properties of the spectral decompositions of the K-matrix 
(p8|), ( p9| ) as we shall explain now. Projecting (|ll]), jl^ ) 
onto transversal mode A we find 



,(N) 



y k n .X 



d_ 

dx 



*S(0) = v / ^«n- (38) 



(D) 



Since the coefficients c^Ta ^ are according to (| 



of order 

2/VA which is independent of n, we have from (|3lf), (|32] ) 



W. 



(N) 
A 



(D) 



(39) 



For a 2D billiard k n ,k ny \ = 0(y/n) (n — > oo) such that 
the terms in the infinite spectral sum ( |27| ) decay asymp- 
totically as n~ 3 / 2 for Neumann and as n -1 / 2 for Dirich- 
let b.c. Hence the convergence is absolute for Neumann 
b.c. while (^9|) converges only conditionally. As a conse- 
quence, the numerically necessary cut-off in the summa- 
tion over the internal states n introduces large errors for 
Dirichlet b.c. 

For the following considerations we will concentrate on 
Neumann b.c. As already mentioned, the choice of the in- 
ternal and external part of the system is not unique, since 
an arbitrary part of the ideal waveguide can be consid- 
ered part of the internal system. Increasing the length of 
the waveguide included, we decrease in fact the influence 
of the evanescent modes since they are exponentially van- 
ishing inside waveguide. We have checked this relation 
and evaluated the time delay functions also for various 



waveguide parts included into the internal system. The 
results remain practically unchanged regardless on the 
length of the included part. This demonstrates the small 
influence of the evanescent modes on the resulting Hei- 
delberg ^-matrix. 

Knowing the coupling matrix W and using the rela- 
tion Eq. (|3|) we evaluated numerically the value of the 
coupling constant g. The obtained result is in an ex- 
cellent agreement with the estimated value (|36| ) for all 
considered types of the waveguide attachment leading to 
g « 0.98. We have evaluated the coupling constant g also 
for a different shape of the billiard [ fl6| obtaining similar 
values for g. It has to be stressed that the estimate ( |36| ) 
was obtained using semiclassical arguments. Our calcu- 
lation shows, however, that it remains valid even in the 
deep quantum region. 

The eigenvalues of H e g (Eq. ||) are usually interpreted 
as the resonance poles and are used for the study of 
the statistical properties of resonances in open quantum 
chaotic systems . To check the validity of this approach 
we have evaluated the resonance poles of the system in- 
dependently using the complex scaling method [M that 
provides a direct access to the positions of the poles of 
the analytically continued S'-matrix. The obtained re- 
sults were compared with the eigenvalues of the effective 
Hamiltonian (|^). However, a direct comparison is ob- 
scured by the fact that the coupling matrix W is in fact 
energy dependent. In the standard random matrix ap- 
proach the coupling matrix is treated as being energy 
independent — a simplification that is well justified inside 
a small energy interval. To mimic this situation and to 
minimize the influence of the energy dependence of the 
coupling matrix W we have compared the eigenvalues of 
.ffeff with the directly evaluated resonance poles always 
within a small energy interval the center of which was 
equal to the energy used to evaluate the coupling matrix 
W. 
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FIG. 4: The complex eigenvalues of H c s (+), zeros of I + 
i7f( N ) ( x ) anc j resonance poles obtained by the complex scaling 
method (o). 
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The energy dependence of W can be taken into ac- 
count more precisely using the relation (|l|) and evaluating 
the Heidelberg resonance poles as zeros of the function 
/ + iK^'(E). We have evaluated the Heidelberg reso- 
nance poles using both of the above described approaches 
and compared the results to the resonances obtained by 
complex scaling. The results are shown in Fig. ^. From 
this figure we see that for narrow resonances the eigen- 
values of effective Hamiltonian H c g represent a good ap- 
proximation to the resonance poles of the system and 
the energy dependence of the coupling matrix W can be 
omitted. For broad resonances the situation changes and 
the complex eigenvalues of the effective Hamiltonian have 
nothing in common with the directly evaluated resonance 
poles. This discrepancy can be explained as follows: A 
resonance localized at Eji = E + iT represents in fact 
a collective mode of all bound states E n of the internal 
Hamiltonian Hi n that are located inside the energy in- 
terval « (E — r, E + T) . For a broad resonance with T 
significantly larger then the mean spacing between the 
bound states E n the number of the internal states to be 
included into H c g must be very high. 

V. CONCLUSIONS 

To summarize, we have shown that the Heidelberg ap- 
proach to scattering, which is the basis of many impor- 
tant random-matrix results on quantum chaotic scatter- 
ing, leads to reasonably good agreement when compared 
with the results of a direct calculation of the S-matrix in 
the case of billiards with Neumann boundary conditions. 
We have explained the somewhat unexpected finding that 
the accuracy is much worse for Dirichlet boundary con- 
ditions, while it does not depend very much on other 
possibilities of varying the auxiliary closed system used, 
such as the position of the attached waveguide. 

Even for Neumann boundary conditions the effec- 



tive Hamiltonian H c g based on the Heidelberg approach 
seems to be not very well suited for the computation of 
broad resonances of the system. This result does not con- 
tradict the fact, that also for systems with time-reversal 
symmetry the statistical properties of billiard resonances 
follow the predictions of random-matrix theory based 
upon the effective Hamiltonian approach quite well [ p3[ , 
because our test goes way beyond a purely statistical 
analysis. 

Moreover, we have shown numerically and with semi- 
classical arguments that the mean dimensionless cou- 
pling for a chaotic billiard is the critical value g = 1 — 
irrespective of the precise form of the billiard, the size of 
the attached waveguide and other details of the model. 
Interestingly, in our model an effective coupling near g = 
1 is observed already deep in the quantum regime. Never- 
theless, fluctuations around the mean value g = 1 should 
in general be largest for small energies and can possibly 
result locally in overcritical coupling. This might be an 
explanation for the observed resonance trapping in bil- 
liards HH0||. 

Our results concerning the value of the effective 
coupling for chaotic systems are not restricted to bil- 
liards and apply, e. g., to quantum graphs as well. In 
these systems, a systematic way to achieve overcritical 
coupling for many states is to modulate the density of 
states, e. g. by considering systems with band spectra 
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